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Abstract 

We study a class of nonlinear Hamiltonians, with applications in quantum optics. The 
interaction terms of these Hamiltonians are generated by taking a linear combination of powers 
of a simple 'beam splitter' Hamiltonian. The entanglement properties of the eigenstates are 
studied. Finally, we show how to use this class of Hamiltonians to perform special tasks such as 
conditional state swapping, which can be used to generate optical cat states and to sort photons, 
keywords: quantum optics, nonlinear optics, quantum information 

1 Introduction 

The consideration of nonlinear optical processes within quantum optics has led to the study 
of many important physical phenomena. These nonlinear optical processes have found many 
diverse applications, such as generating entangled photons using parametric down conversion 
[TJ or creating optical routers using intensity dependent properties of an optical fibre [H [3] . The 
growing field of quantum information has added to the interest in nonlinear optics. This is due 
to the fact that nonlinear optics can be used to perform certain operations that are not possible 
using only linear optics. For example, a class of nonlinear optical Hamiltonians, which have 
been used to model four wave mixers [H 03 [HI [7] , have been shown to allow one to prepare a 
macroscopically distinguishable superposition of quantum states. 

In this paper we will consider a class of nonlinear Hamiltonians. The Hamiltonians will be 
functions (e.g. a polynomial) of a quadratic Hamiltonian, Hq, that leads to linear differential 
equations in the field mode operators. Hamiltonians of this type have been used to describe 
optical beam splitters. The nonlinear Hamiltonians will thus commute with Hq. Hence the 
task of obtaining the eigenvectors and eigenvalues of the nonlinear Hamiltonians will reduce 
to diagonalising the simple quadratic Hamiltonian. It will be shown that the eigenvalues and 
eigenvectors are related to a class of orthogonal Krawtchouk polynomials [8] . We thus establish 
that this particular class of nonlinear Hamiltonians are exactly solvable. The connection between 
orthogonal polynomials and nonlinear optical process has been studied previously and general 
results have been obtained [HI QUI EE] • In particular, the mathematical results that we present 
in section [5] is a special case of the general theory presented in [5] . In [9] the general method 
is described for solving two boson systems in quantum optics via of orthogonal polynomials 
systems. It is possible to study many of these problems using different mathematical techniques. 
In particular, beam splitters have bean studied extensively p~21 E2 [14] . The benefit of using 
orthogonal polynomials, however, is that there exists extensive literature on their properties, 
which one can exploited for any given problem, see for example [51 ITS]. 

The organization of the paper will be as follows. In section[2]we will diagonalise the quadratic 
Hamiltonian. This will be achieved by first showing the equivalence between diagonalising the 
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Hamiltonian and solving a particular recurrence relation. The solutions of the recurrence relation 
are the Krawtchouk polynomials. The properties of these polynomials will be used to obtain 
the spectrum of the Hamiltonian. In section 3 we study the entanglement properties of the 
energy eigenstates. The properties of the nonlinear Hamiltonians will be investigated in section 
[4] In particular we will show how the Hamiltonians can be used to perform a conditionally swap 
operation on the two modes, when the modes are both prepared with all the photons initially 
in one mode. We will also show how this enables one to prepare optical Schrodinger cat states. 
The dynamics of the nonlinear Hamiltonians is studied in further detail in section Finally, 
we discuss our results in section [5J 

2 Diagonalising the quadratic Hamlitonian 

In this section we apply the method investigated in [9]. Suppose we have an electromagnetic 
field that has two modes 1 and 2. Associated with these modes are the creation and annihilation 
operators and aj, where j = 1,2. These operators obey the standard canonical commutator 

relations [aj,a^] = i and [dj,a,k] — 0, for j ^ k. Suppose further that the two modes interact. 
If this interaction is quadratic, then one possible form for the interaction Hamiltonian is 

H = 7a[a 2 +7*o|ai. (1) 

where the free evolution terms are ignorecfj] Hamiltonians of the form {l} are be used to describe 
beam splitters. A simple class of nonlinear interaction Hamiltonians can be constructed by 
introducing some nonlinear function f{x) and forming a new Hamiltonian H = f(H ), where 
Hq is defined in equation ([T]). A simple example of this would be f(x) a polynomial e.g. 

H = ^H Q + \{H a f. (2) 

The nonlinear Hamiltonian described in equation ^ has been studied before in connection with 
various nonlinear optical phenomena, see for example [H [5] . 

The quadratic Hamiltonian Hq conserves the total number of photons, i.e. 

[n,H o ] = 0, (3) 

where h — h\ + n 2 . The class of nonlinear Hamiltonians H = /(Hq), where f(x) is some polyno- 
mial, will also conserve the total number of photons. By construction we see that [Hq, H] = 0. 
We can thus obtain the eigenvectors of H by diagonalising Hq. The fact that both Hamiltonians 
conserve the total photon number suggests that we should diagonalise the Hamiltonian (fTJ) in 
subspaces with fixed numbers of photons. If we have a subspace where the total number of 
photons is constant and equals M , then the dimensions of this subspace will be M + 1 and any 
state belonging to the subspace can be written in the form 

M 

|<A M }i2 = ]T£n|M-n,n) 12 , (4) 

71=0 

where \M — n, n)i 2 represents the Fock state \M — n}i (g) \ri)2- Let \E M )i 2 denote an eigenstate 
of (fT]), hence H \E M )i 2 = E\E M )i 2 . For the sake of convenience we shall express the eigenstate 
in the following form 

AI 

\E M ) 12 = c^e~ m9 \M - n 7 n} 12 , (5) 

where g = arg(7), i.e. 7 = |7|e Ifl . Acting on equation |5]) with Ha yields the expression 
H \E M ) 12 ^ 1 Y J ^Jn{M -n+ T)e- in *c*f\M -n + l,n- 1) 12 

n 

+7* J2 V(n + l)(M~n)e' m9 c^\M -n,n+ 1) 12 = E\E M ) 12 . (6) 



1 One could imagine that we are working in an interaction picture. 
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By re- arranging this expression one can show that the coefficients c^f , obey the following three 
term recurrence relation 



E 



— cf = v /( n + l)(M-n)cf +1 + \/n(M — rc + 1) 
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(7) 



The recurrence relation can be solved to find both the eigenvalues and the coefficients that 
appear in the Fock basis expansion of the energy eigenvectors. The mathematical properties 
of the recurrence relation ([7J enable us to determine many properties of the eigenvectors and 
eigenvalues of energy. For example, if equation ([5]) is an eigenvector of Hq corresponding to 
the eigenvalue E, then so is the vector | — E M )i 2 = J2 n ( — 1)™ ex P( — ^ n 9) c n \M — n, n)i 2 , which 
corresponds to the eigenvalue —E. The proof of this fact follows simply by replacing the terms 
Cn i m 0i with (— l) n c„ . We thus see that for every positive eigenvalue of energy there must 
exist a corresponding negative energy eigenvalue. 

The coefficients c^ 1 will be functions of E. If we formally solve the recurrence relation (J7J, 
then we find that c^f (E) is a polynomial of degree n, in the variable E. The fact that the total 
number of photons is M means that cj^ I+1 (E) = 0, which means that the roots of the polynomial 
cjf I+1 (E) correspond to the eigenvalues of (JTJ) . To help explain the last fact consider the matrix 
representation of equation |7} 
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(8) 



M+i is the term that we would find in front of c\\ +l in the recurrence relation ([!]). 



where 

X (E) = 0, then the column vector c M will be an eigenvector of the 
Hamiltonian, corresponding to the eigenvalue E/\j\. An important point to note about equation 

and c%Li are positive. This fact enables us to 
For a recurrence relation of the form 



J7]) is that the factors that multiply both c^+i 
make use of a theorem of Favard's [17] 



xP„(x) = b n+1 P n+1 (x) + b n -iP n -i(x), 



(9) 



where the sequence {b„} is non-negative, there exists a linear functional C, which has the 
following properties 



£[1] = 1, 

C[P m {x)P n {x)\ = Oform^n, 
C[P*(x)] f 0. 



(10) 



Applying this theorem to ([7]) leads to the conclusion that the set of polynomials {^(x)}, will 
be orthogonal to each other with respect to the linear functional C. A standard property of 
orthogonal polynomials is that the roots of an orthogonal polynomial are distinct |16j . This 
last result together with the fact that the eigenvalues of Hq are the roots of the polynomial 
Cm +1 (E) = 0, implies that the eigenvalues of Hq are nondegenerate, when we are confined to 
the M photon subspac^E There is an extensive literature on orthogonal polynomials and their 
properties [HEStlB]. One important example are the Krawtchouk polynomials, which obey the 
following three term recurrence relation [8] 



xK n (x; p, M) = ^/p(l-p)(n + l){M -n)K n+1 (x; p, M) + {p[M - 
+VP(1 - p)n(M ~n+ l)K n ^{x- p, M), 



■n[l-p])K n (x;p,M) + 
(11) 



where < p < 1. If we set p = 1/2, 2x - M = E/\j\ and K n (x; 1/2, M) = c*f(E), then 
equation pip reduce to equation ([7]). The coefficients c^f (E) are thus Krawtchouk polynomials. 



2 The spectrum of Ho on the whole Hilbert space will have degeneracies due to the fact that the spectrum of Hq 
on the different subspace is not disjointed and thus can share eigenvalues. 
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The form of the Krawtchouk polynomials, for this choice of parameters, is given by 



A-„(,; 1/2. M) . (-!)•*„(„ 1/2, M)^^L_ g (12) 

where (j/)^ = y(y + l)(y + 2) . . . (y + k — 1). The roots of the Krawtchouk polynomials can 
be determined using (fT2)) . It can easily be shown that the roots of Km+i(x; 1/2, M) = are 
x = 0, 1, 2, M, which implies that the energy eigenvalues are 

E = -M| 7 |, (2 - M)| 7 |, (4 - M)| 7 |, M| 7 | (13) 

By diagonalising the Hamiltonian ([I]) we have solved the system's dynamics and can analyze 
relevant quantum effects. 

3 Photon distribution and entanglement properties of the 
energy eigenstates 

In the previous section we found the analytic solution of equation ([7]), which in turn gives us 
the form of the eigenvectors of the Hamiltonian ([T]). Nevertheless, it is informative to give some 
simple examples of the eigenvectors of ([T]). For the simplest case, when M = 1, the eigenvalues 
are just E — ±| 7 | and the eigenvectors take the form \E±) = (|10) ± exp(— ig)\01))/^/2. When 
M = 2, we find that E = ±2| 7 |, and 

|£ = ±2| 7 |) = i(|20)±V2 e ^|ll>+ e - 2 ^|02)), 

\E = 0) = i=(|20)- e -^|02)). (14) 

As M increases the form of the eigenvectors will generally become more complicated. One 
situation, however, where the eigenvectors have a simple form is when E = M| 7 |, i.e for the 
eigenvectors corresponding to the maximum eigenvalue of Hq . In this case we find 



k=0 



( 1V , M Ml ^ n\ k _ M Ml 

~ [ } C ° y n \(M- n y.^ k\(n-k)\ [ ' ~ \]n\(M-n)V [ ' 

If Cff is chosen to be 2~ M / 2 , then ^2 n \cff (M\^\)\ 2 = 1. A consequence of this result is that the 
photon probability distribution, |c^ f | 2 , will be a binomial distribution when E = M| 7 |. In the 
limit of large M, |e^(M| 7 |)| 2 can be approximated by a Gaussian. 

The probability distribution |c* f | 2 , of finding n photons in the second mode, depends only 
on \E\. By this we mean that \cff (E)\ 2 = \cff(—E)\ 2 . This follows directly from the result 
in section [2] that c^f(—E) = (—l) n c^f(E), which can be verified using the recurrence relation 
J7]). The behaviour of the probability distributions \c^f(E)\ 2 , for \E\ < M| 7 |, can be studied 
by numerically evaluating (fT2| . Figure Q] shows | 2 plotted with M = 38, for the four largest 
eigenvalues of Hq. We see that each time E is decreased by 2| 7 |, an additional peak appears 
in the probability distribution. Numerical investigations show that this behaviour is a generic 
feature and occurs for each value of M . Figure [1] shows that the probability distribution has the 
symmetry, \c^\ 2 = |c*^_„| 2 . This means that the probability of finding n photons in the first 
mode is the same as finding n photons in the second mode. This result is physically reasonable 
due to the symmetric nature of the Hamiltonian ([T]), with regards to the two modes. 

If we calculate the explicit forms for the energy eigenvectors, then we see that the two modes 
are entangled, i.e they cannot be express in a factorised form. Entanglement is one of the defining 
characteristics of quantum mechanics and has great importance in the foundations of quantum 
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Figure 1: Four different piots of the probability distribution for finding n photons in the second 
mode. Each plot is for the situation where M = 38 and I7I = 1, however in each successive plot 
the energy is decreased by 2\y\. 
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Figure 2: A plot of the entanglement, S en t, against the number of photons M, for different values 
of energy E, with 7 = 1. 



mechanics and in the field quantum information [18] . For this reason it is of interest to quantify 
the entanglement present in the energy eigenstates. A common approach to quantifying the 
entanglement in a bipartite pure state is to calculate the von Neumann entropy in the reduced 
state [TBI HH1 HO] ■ The von Neumann entropy of a density operator p is defined to be 

S[p]=-Tr[plog(p)}, (16) 

where the logarithms are taken to the base 2. The reduced states of the bipartite state |V'}i2 
are define to be pi = Tr 2 [|i/')('0|] and f>2 = Tr i[\ip) (tp\], where Tr^ denotes the partial trace over 
the subsystem j. We thus take use the following quantity to measure the entanglement [19] 

S ent {\il>)y>\)=2S(p 1 ) = 2S(fi 3 ). (17) 

S'entdV'KV'l) — when our bipartite state is separable. In addition to this, we can see that 
equation (| 1 7[) assumes its maximum value when our bipartite state is a maximally entangled 
state, i.e. the reduced states are proportional to the identity operator. If we express the energy 
eigenvectors in the form given in equation ([5]), then the eigenstates are automatically in the 
Schmidt form [18], where the basis states {\M — n, 71)12} now act as our Schmidt basis. From 
this observation is is straightforward to show that Sent equals twice the entropy of the probability 
distribution Ic* 1 ^)! 2 , i.e 

S ent = 25(|cf | 2 ) = -2^|c^| 2 log(|cf | 2 ), (18) 

n 

where the logarithms are taken to the base 2. We can now evaluate S en t as a function of the total 
number of photons M. Figure [2] shows S en t plotted against M, with the energy E either fixed 
or set equal to the largest possible value given the choice of M. We see that as M increases so 
does S en t- The fact that the entanglement increases with M is not surprising as the dimensions 
of our systems Hilbert space are M + 1. Increasing the number of photons thus gives us more 
terms in the Schmidt decomposition of the eigenstates \E M )i2- We can also see that the zero 
eigenstate, i.e. the state \E M )i2, with E = 0, is less entangled than the states with E = 1 or 
E = 2. By examining equation ([7]) is can be seen that for E — 0, we would require that every 
term equal zero for n odd. The zero eigenvector thus has roughly half the number of Schmidt 
terms as the other eigenvectors do. Consequently, the zero eigenvector is not as entangled as 
the eigenvectors that correspond to E 7^ 0. 

It is interesting to see how the entanglement changes as we vary E, with the total number of 
photons now fixed. Figure [3] shows S en t plotted against the eigen-energy, E, for several different 
values of M. It can be seen that when E is large, the entropy decreases. The fact that there 
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Figure 3: Entanglement, S en t, plotted against the energy E, for different values of the photon 
number M. 

is a general trend for the entanglement to decrease as E increases can be explained by looking 
at how the probability distribution, \cff(E)\ 2 , changes when E varies. Figure [1] shows that 
as we decrease E additionally peaks appear in the probability distribution. Fhis means that 
the probability becomes more spread out and thus the entropy increases. Thinking in terms 
of correlations we see that for E equal to M|7|, we have a single maximum in the probability 
distribution |c^| 2 . The Schmidt basis states \M — n,n) 12 will only have coefficients with large 
values for the states with n close to M/2. When E is decreased we obtain several maxima in 
the probability distribution \cff (E)\ 2 and thus we have more Schmidt basis states with larger 
coefficients and hence the entanglement increases. 

An interesting feature of figure [3] is that there are points where the entanglement increase 
slightly as E increases. One obvious example of this is when M is even and E increases from 
zero to 2 (7 = 1 in figure [3]). In this case the increase is due to the anomalous nature of the 
zero energy eigenvector, i.e it has roughly half the number of nonzero Schmidt coefficients as 
the other energy eigenvectors. The other instances when the entanglement increase occur due 
to the fact that P^f '(E) can sometimes broaden when E increases. Sometimes this broadening 
will compensate for the loss of peaks in P r A4 (E) and thus the probability distribution is more 
spread out, which leads to a slight increase in the entropy. 

4 Nonlinear Hamiltonians and conditional state swapping 

In section[2]we obtained the eigenvalues and eigenvectors corresponding to the quadratic Hamil- 
tonian ([!]). These results can now be used to determine the properties of a certain class of 
nonlinear Hamiltonians, which are constructed from Hq [an example of this is given in equation 



Let f(y) be a polynomial in the variable y, it is clear that [Hq, /(Hq)] = 0, and thus 
f(H )\E M )i2 = f(E)\E M ) 12 . Similarly, the operator h = h\ + h 2 , obeys [h,H ] = and thus 
f(h)\E AI ) 12 = f(M)\E AI ) i2- From this we see that the general Hamiltonian 




(19) 



k 



can be diagonalised using the eigenvectors {|-E M )i2j-. We should note that the spectrum of 
equation (Tl9|) will generally be degenerate and thus we can construct several different sets 
of eigenvectors for the Hamiltonian. One important practical consideration for our nonlinear 



Hamiltonians is that the magnitude of the terms a k should decrease as k increases. In practice, 
if we want the higher order terms in our Hamiltonian to give an non-negligible contribution to 
the evolution of our system, we would have to increase the intensity of our light beams. This 
would correspond to M being large. Even when we increase the intensity of the optical beams 
that enter the nonlinear media, it will still be difficult to excite higher order terms. For this 
reason we shall look only at examples were H is at most quadratic in H , i.e. Hamiltonian with 
a form given in equation @. 

Using our previous results we can easily determine the time evolution of a system described 
by a Hamiltonian of the form given in equation (|19[) . If at time t — the system is prepared in 
the state 1^0)12, then at a later time t, it will be in the state \ipt)i2 = exp(— iHt)\t()o)i2, where 
h = 1. The problem of determining the evolution thus reduces to expressing the state \ipo)i 2 
in terms of the basis {\E M )i 2 }. In quantum optics, states are often expressed in terms of the 
Fock basis. It is thus essential to be able to express a given Fock state in terms of the energy 
eigenvectors. This can be achieved using the completeness relation, ^2 k \E k )(E k \ = 1, which 
enables us to obtain that 

\M-n,n) 12 =J2i2( E k\M-n,n) 12 \E k ) 12 = e in9 ^c™(E k )\E k ) 12 . (20) 

k k 

We now give examples of how the preceding theory can be used to construct Hamiltonians 
that perform specific tasks. In order to simplify the following calculations we set 7 real. Suppose 
we have the Fock state \M, 0)12, which we want to transform to the state |0, M)\ 2 , i.e. we swap 
the two modes. From equation (|20|1 it is clear that if we want to perform this task, then we 
must change the coefficient Cq 1 (E k ) to c^(Ek). Using equation (TT2"]) together with some simple 
algebra leads to the result that 

c%(E k ) = (-l) M+ *><g(E k ), (21) 

where x k — (E k + M|7|)/(2|7|). This result means that when M + x k is odd, cj^(E k ) = 
—C^'(Ek), while for M + x k even Cf^(E k ) = Cq 1 (E k ). Suppose that we wish to construct 
a Hamiltonian that will peform the swap operation \M, 0)i 2 — > |0, M)\ 2 for all values of M. 
There is an infinity of different Hamiltonians, which achieve this task. An example of a simple 
Hamiltonian that performs the swap operation is 

^=^(3M" + ffo), (22) 

where exp(— iHt)\M, 0)12 = |0, M)\ 2 . This Hamiltonian is linear in Hq, we could, however, 
construct alternative nonlinear Hamiltonians that perform the swap operation. An interest- 
ing point to note is that if we had some state \tp) = X)n£™l n )' then the Hamiltonian de- 
scribed in equation (j22|) can be used to enact the transformation |^)i|0)2 — > |0) 1 1 -0} 2 , i.e. 
exp(-iffT)|^)i|0) 2 = |0)i|V>>2- 

In the previous example the swap operation was insensitive to the total number of photons. 
It is, however, possible to design Hamiltonians that will perform a swap operation only if the 
number of photons is either even or odd. Suppose now that we want the swap operation to occur 
only if the total number of photons is even, i.e |M, 0)12 — > |0, M) 12 only when M is even. Let 
e^ 1 denote the energy eigenvalue of our nonlinear Hamiltonian, corresponding to the eigenvector 
\E^)i 2 . In order to perform the swap at time t — r we require that exp(— ie^fr) = (— 1) M+X , 
but only when M is even. A Hamiltonian that achieves this is 

n ( , I7L2 , 1 £j2 



One can easily verify that 



H = W\ m+ 4* + W\ Ho] - (23) 



- \ 1 , f 10, M) 12, for M even, , , 

exp-^rM,0 12 = ,2 ' ' (24) 

i\M, 0)i2, for M odd. 
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The Hamiltonian (|23[) thus has the property that for a state \ijj), the transformation |0)a — * 
1 0)i 1^)2 can only be realized if can be expressed as a superposition of Fock states consisting 
of an even number of photons. 

We have described how the Hamiltonian given in equation (|23p . can be used to perform a 
conditional swap operation. The Hamiltonian can also be used to perform other useful tasks, 
such as preparing two optical cat states form a single coherent state. An optical cat is a state 
of the form K(\a) + | — a)), where K is a normalisation constant and | ± a) are coherent states 
[?]. The significance of these states can be seen by recalling that when \a\ is large, coherent 
state behave like a classical monochromatic optical fields. Cat states are thus a superposition 
of two classical object^. The Fock state representation of a coherent state is [22] 

where a is a complex number. If we take our initial state to be |a)i|0)2, then using the 
Hamiltonian (|23l) will yield the following evolution 

ex P (-^)|a) 1 |0) 2 = |0>x ( Ia)2+ j~ a)2 ) |0), (26) 

We have thus obtained a superposition of two optical cat states. Theoretical schemes for obtain- 
ing a superposition of two optical cat states have been proposed previously [H [5] [7] . In each of 
the earlier schemes, the Hamiltonians that were used were different from that given in equation 

HMD- 

As we have shown, it is possible to construct Hamiltonians that perform a swap operation 
conditioned on whether the number of photons is even. It is also interesting to consider a swap 
operation conditioned on other factors. For instance, we may want to not swap the state if it 
has a particular number of photons, while allowing the swap to occur if the number of photons 
differs by one. We now give the form of a Hamiltonian that performs this task. As before, let 

denote the energy of the nonlinear Hamiltonian. We require that exp(— ietfr) = (— \) M+X 
if M = N ± 1, but that exp(-ie* / r) ^ (-1)^+* for M = N. It can easily be verified that the 
following Hamiltonian satisfies the stated conditions 

H=^^{3\-f\h 2 + hH -3N\l\n~NHo). (27) 

This leads to the following dynamics, exp(— iHr)\N, 0)12 = |iV, 0) 12, while exp(— iHr)\N ± 
1, 0)12 = |0, iV± 1)12. It is straightforward to modify the Hamiltonian given in equation ([27]) for 
situations where we want to swap states other than |iV± 1, 0)12. An example of this is when we 
want to swap states such as \N ± 2, 0)12. To achieve this we can simply multiply equation ([27]) 
by a factor of 1/2, i.e. we take our Hamiltonian to be H/2. We shall now show how one can use 
the Hamiltonian (|27p together with our previous nonlinear Hamltonians, so as to discriminate 
between photon number states. 

Suppose we have four modes that are initially prepared in the state 

|M)i|0) 2 |0> 3 |0>4, (28) 

where M < 4. An important problem to address is whether we can determining the value of M. 
The majority of currently available photon detectors are not capable of reliably discriminating 
between different photon number states. For this reason it would be interesting if we could act 
on the state (|28|) , so that the final position of the mode that is not in the vacuum is dependent 
on the number of photons. In particular, we could arrange for the state to evolve such that the 
M-th mode contains M photons, while the other modes are in the vacuum state. We shall now 
outline a scheme that enables one to performs this task. 

3 The name 'cat state' is a reference to the thought experiment developed by Schrodinger [23] . 
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At time t = our system is prepared in the state Q28p . We arrange for modes 1 and 2 to 
interact in a manner governed by the Hamiltonian (|23[) . This means that at time r the photons 
will have been transferred to the second modes if M is even, or remain in the first mode for M 
odd. We now arrange for the first and third modes to be coupled via a Hamiltonian of the form 
H/2, where H is given by equation (|27|) and where the N in equation (f27|) equals 1. This will 
result in the first and third mode being swapped when we have three photon. At the same time 
we shall also couple the second and fourth modes using a Hamiltonian of the form H/2, where 
H is given by equation (|27p . with N = 2. This will result in the second and fourth modes being 
swapped if the number of photons is four. The effect of coupling the modes in the manner that 
we have described is that our system will be left in the state 

i|l)i|0)a|0) 3 |0) 4 , for M=l, 
|0)i|2) 2 |0) 3 |0) 4 , for M=2, 
|0)i|0) 2 |3) 3 |0) 4 , for M=3, 

|0)i|0) 2 |0) 3 |4) 4 , for M = 4. (29) 

We could now put photon detectors at each mode so as to determine which mode contains 
photons. This in turn would allow us to determine the number of photons that we have. 
Alternatively we could perform further operations on the individual modes. In addition to 
this, the above scheme can be used to prepare multiphoton entangled states, provided we can 
prepare a superposition of a finite number of Fock states, in a single modes. For example, 
we could transform the separable state (a\l) + 6|2))|0) into the entangled state a|10) + 6|02). 
The various procedures that we have outlined could find applications in quantum information 
processing and quantum control. 

5 Dynamics of the nonlinear Hamiltonians 

In the previous section we constructed nonlinear Hamiltonians that allowed use to perform 
certain types of operations on the two coupled modes. In this section we will investigate further 
the dynamics induced by the three Hamiltonians given by (p2| . (j23f and (j27|) . A natural quantity 
to study is the number of photons in each mode. For this reason we shall investigate how the 
expectation value for the number of photons in the second mode changes, i.e. (n 2 ). The initial 
state of our system will be taken to be |M, 0)i 2 , and thus (n 2 ) = at time zero. 

We shall now study the dynamics associated with the quadratic Hamiltonian (|22[) . In figure 
[4] we see the expectation value (n 2 ) plotted as a function of time, for the initial states 1 10,0) 12 
and 1 1 1 , 0) 12 , with 7 = r = 1. The expectation value oscillates with period of two (i.e. 2r), 
the period is thus independent of the number of photons. Another interesting feature of figure 
[4] is that at time t = 1, (fi 2 ) = M, which is consistent with the state having been transformed 
to |0, M)i2 at that time. The oscillatory dynamics shown in figured] is also observed for other 
values of M, and does not depend on whether M is even or odd. One final interesting feature 
of the Hamiltonian ([22]) is that we obtain oscillations for different choices of the initial state. 
For example, if the state was initially \M — n,n), then we would still observe oscillations of 
(n 2 ), with period 2r. The minimum and maximum values of (n 2 ) would, however, respectively 
change to (n 2 ) = n at t = and (n 2 ) = M — n at t = r. Such behaviour is consistent with the 
Hamiltonian (|2"2"]l being able to swap the two modes. 

The Hamiltonian (|23[) was designed to perform the transformation described in equation 
([24]) . This selective transformation has similarities to the behaviour of a two mode nonlinear 
directional coupler [21 [3] . This is an optical device which consists of two optical fibers that are 
coupled in such a way that the dynamics depends on the intensity of the light |3J . One obvious 
difference between these two nonlinear systems is that the conditional dynamics associated with 
the Hamiltonian |23]) . dependent on whether the number of photons is even or odd and not on 
how many photons there are. Figure [5] shows a plot of the expectation value (n. 2 ), against the 
time, for the initial states 1 10, 0) 12 and 1 11, 0) 12- We see that the dynamics is very different for 
the two different initial states. In particular, (n 2 ) never exceeds the value 5.5, for the case when 
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Figure 4: A plot of the expectation value (^2}, against time, with 7 = r = 1 and where the initial 
state is \M, 0}n- 

M = 11. This is in contrast to the situation for M = 10, where we find that (77,2) = 10 at t = 1 
(i.e. t = t). In both cases we see that the behaviour of (712) is periodic, with period t = 2t = 2. 
When M = 10, one can observe that the expectation value (^2) temporarily plateaus at the 
value of 5, which is half the total number of photons. 

It is interesting to investigate the dynamics of the Hamiltonian (|24|) for initial states, which 
have photons in both modes, i.e \M — n, n)i2, with n^O. Figure (J6j) show (^2) plotted against 
time, for the initial states |S, 2) 12 and |9, 2) 12 - In both cases the minimum of (71,2) is now 2. 
For the case when M — 10, we see that at t = 1, (112) = 8, which is consistent with the state 
undergoing the transformation |8, 2)12 — > |2, 8) 12 - When the initial state is |9, 2) 12 and thus 
M = 11, {712) is always less than 9, and thus there is never a point where we transfer the 9 
photons that were initially in the first mode, to the second mode. 

The Hamiltonian ([27]) was designed so that it would perform the transformation \N ± 
l,0)i2 — > |0, -/V ± 1)12, while leaving a state with exactly TV photons unchanged. It is clear 
that if our system is initially prepared in the state | TV, 0)12, then the system will stay in that 
state. In fact, any state with exactly TV photons will remain unchanged under the action of (|27|) . 
If our system has TV ± 1 photons then the dynamics of the Hamiltonian (|27|) will be identical to 
that of (22|). 

6 Conclusions 

We have investigated a class of nonlinear Hamiltonians. These Hamiltonians were constructed 
from a simple beam splitter Hamiltonian, equation {T]), and from the photon number operator 
n. Due to the pivotal role played by the beam splitter Hamiltonian, |T]), we investigated its 
properties thoroughly. The eigenvalues and eigenvectors of equation |T]) were obtained using 
the theory of orthogonal polynomials. In particular, it was shown that the eigenvalue problem 
associated with equation (TT]) was equivalent to solving the three term recurrence relation as- 
sociated with a class of orthogonal polynomials known as Krawtchouk polynomials. We were 
thus able to show that the dynamics of the beam splitter Hamiltonian and the class of nonlinear 
Hamiltonians, are exactly solvable using Krawtchouk polynomials. 

Nonlinear Hamiltonians were then constructed from the beam splitter Hamiltonians and the 
total photon number operator. The dynamics of these Hamiltonians enabled novel tasks such as 
state swapping and conditional state swapping, to be performed. The formalism that we have 
introduced and the Hamiltonians that we have studied, have applications in quantum optics 
and quantum information. In particular, they can be used in the areas of state preparation and 
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Figure 5: A plot of the expectation value (^2), against time, with 7 = r = 1 and where the initial 
state is \M, 0)12. It can be seen that for M = 11 the Hamiltonian acts to trap some of the energy 
in the first mode, while for M = 10 the energy is fully swapped between the two modes. 




Figure 6: A plot of the expectation value (^2), against time, with 7 = r = 1. In plot (a) the initial 
state is 1 8, 2)12, while in plot (b) the initial state is 1 9, 2) 12 . 
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in quantum control. For example, we have shown how the nonlinear Hamiltonian (|23p can be 
used to prepare a superposition of Schrodinger cats states for a coherent state. In addition to 
this we have shown how we can control several modes so as to sort photons into one particular 
mode, where the mode that there are transferred to depends on the total number of photons. 
This was achieved by turning on couplings between two different modes, where the interactions 
were described using the nonlinear Hamiltonians (f2"3")l and (f2"T)) . While we have discussed this 
work in terms of optics and photons, the results would apply to any system of bosonic particles. 
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